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SUMMARY 

A procedure was developed and is described which can be used to computa- 
tionally simulate the cyclic behavior of high temperature metal matrix compo- 
sites (HTMMC) and its degradation effects on the structural response. This 
procedure consists of HTMMC mechanics coupled with a mul ti factor-i nteraction 
constituent material relationship and with an incremental iterative nonlinear 
analysis. The procedure is implemented in a computer code which can be used 
to computationally simulate the thermomechanical behavior of HTMMC starting 
from the fabrication process and proceeding through thermomechanical cycling, 
accounting for the interface/interphase region. Results show that combined 
thermal /mechani cal cycling, the interphase and in situ matrix properties have 
significant effects on the structural integrity of HTMMC. 


INTRODUCTION 

Structural components made from high temperature metal matrix composites 
(HTMMC) will be subjected to both monotonic and cyclic thermal and mechanical 
loads. The cyclic loads degrade the composite structural integrity in terms 
of stiffness and strength. What is of Interest is either (1) the degradation 
rate or (2) the value to which characteristic properties have degraded after a 
specified number of cycles. Experimental determination of these degradations 
is frequently prohibitive in cost and may not be tractable because of the 
numerous contributing factors such as constituent material properties, fabri- 
cation process variables, and the continuously changing nonlinear behavior with 
each cycle. 

Recent research by the authors indicates that the degradation can be 
described at the constituent materials level, including the interphase, by 
using a multifactor-interaction relationship of product series form. The cor- 
responding degradation at the composite level is described by using composite 
micro- and macro-mechanics incorporated into a computer code. This computer 
code is used to computationally simulate the cyclic behavior of HTMMC. The 
objective of the present paper is to describe the computational procedure and 
illustrate its application by presenting typical results. 


GLOBAL STRUCTURAL RESPONSE 

A brief description of global structural analysis provides the background 
for the need to computationally simulate the thermal /mechani cal behavior of 
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HTMMC. The process of "computational simulation" is summarized in the follow- 
ing section. A schematic summarizing equations for global structural analysis, 
through which global structural response variables and local variables are 
obtained to assess the adequacy of the design, is depicted in figure 1. Com- 
putational simulation is used to generate the elements (terms) in the struc- 
tural matrices: CM], [C ] , CK] , CD], and {S^}. The vector { F< t ) } Includes the 

loading conditions while {u} includes the calculated global displacements, {a} 
includes the local stresses and u represents the natural frequencies of the 
structure. Variables with the subscript A denote allowable values of the 
respective variables set by the design requirements. Comparable equations are 
used for thermal analysis where composite micro- and macro-mechanics are used 
to generate the structure's thermal properties required for heat transfer and 
for thermal stress analysis. The sections that follow describe the details of 
the procedure to generate the properties needed for local /global thermostruc- 
tural analysis of HTMMC. 


COMPUTATIONAL SIMULATION OF THERMAL/MECHANICAL PROPERTIES 

For the computational simulation to be representative of the physics of 
HTMMC behavior, It must account for the fabrication process starting with the 
constituent materials and ending with the structural component. A schematic 
depicting the fabrication process of HTMMC is shown in figure 2. The integra- 
ted multiscale computational simulation strategy is shown in figure 3. The 
left part of the figure represents the behavior synthesis, while the right 
part represents the response decomposition. In the synthesis the properties 
needed for structural analysis are generated, while in the decomposition the 
local stresses are calculated. 

The dependence of the constituent material properties on time, tempera- 
ture, and stress, is defined in a materials space shown at the bottom center 
of the figure. The ply stresses are used to calculate the local microstresses 
(the stresses In the constituents and in the interface). The simulation (syn- 
thesis and decomposition) is performed using an incremental-iterative nonlinear 
analysis where convergence is enforced at every level (ply, laminate, and 
global) of the simulation. A detailed description of this simulation includ- 
ing requisite equations is found in reference 1. A stand-alone computer code 
has been developed to perform the computational simulations from the constitu- 
ent materials level to the laminate level. This stand-alone code encompasses 
all the capabilities shown within the dashed line in figure 3 and is identified 
as METCAN (Metal Matrix Composite Analyzer). METCAN is patterned after and has 
many of the features of the ICAN computer code (ref. 2). 

The various factors, which influence constituent material behavior and 
ultimately the HTMMC behavior, are incorporated through a multifactor- 
interaction relationship as shown in figure 4 where the rationale for adopting 
this type of relationship is summarized. Additional discussion on this rela- 
tionship and its application to a specific structural component fabricated 
from HTMMC is found in reference 3. The important points to note in figure 4 
are: (1) the individual effects of the thermal and mechanical load cycles 

Ny , N|v) are explicitly included; (2) the temperature T and the stress a 
induced by these cycles are also explicitly included; (3) the combined 
thermomechanical cyclic effects on the thermomechanical behavior of HTMMC' s 
are, therefore, automatically included; and (4) the various effects are 
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normalized with respect to reference properties which are usually taken as room 
temperature properties. 

The important points of the above discussion are (1) METCAN can be used 
to computationally simulate the high temperature nonlinear behavior of HTMMC' s 
from the fabrication process through service and from global structure response 
to local constituent response including the Interface/interphase regions, and 
(2) METCAN can be used as a preprocessor to any global structural (thermal) 
analysis with anisotropic material capability and also as a postprocessor to 
this same global structural analysis. METCAN is a comprehensive and continuo- 
usly evolving computer code. In this respect, it is difficult at present to 
verify all of its predictive capabilities since most of the necessary data do 
not exist. Even so, the code is continuously validated (compared with other 
computational predictions) and verified (compared with available or currently 
acquired experimental data). The room temperature constituent properties used 
to generate the results that follow are given in Table I. Verification compar- 
isons of room temperature stress/strain response for tungsten/copper HTMMC are 
shown in figure 5 for two fiber volume ratios. Similar results for graphite/ 
copper, but only at one fiber volume ratio, are shown in figure 6. The exper- 
imental data was obtained by J. Diaz, Materials Division, NASA Lewis Research 
Center. Note that METCAN predicts higher stress and strain to fracture for 
the graphite/copper which could Indicate fiber damage during the fabrication 
process. Verification and validation comparisons for borsic/ti tanium metal 
matrix composite are shown in figure 7. The agreement is good in all these 
comparisons, except as already noted. Validation comparisons for room temper- 
ature thermal and mechanical properties for a specific graphite fiber (PI 00) / 
copper matrix (Cu) HTMMC with four different fiber volume ratios are summarized 
in table II. These comparisons also show excellent agreement. Though these 
comparisons are admittedly limited, they do indicate that one can computation- 
ally simulate the thermomechanical behavior of HTMMC using METCAN with some 
confidence. 


COMPUTATIONAL SIMULATION OF CYCLIC BEHAVIOR 

The thermal/mechanical cyclic behavior of HTMMC is computationally simu- 
lated by using the following procedure: 

(1) Calculate the microstresses from global thermal and structural analy- 
ses using the respective cyclic loads as inputs to these global analyses. Note 
the global analyses are performed using finite element analysis of the struc- 
tural component (top part (undashed) fig. 3). 

(2) Enter the magnitude of these microstresses and the respective number 
of thermal and mechanical cycles in the multifactor relationship in figure 4. 

(3) Apply the next cycle loads to the structure and/or laminate and simu- 
late the desired behavior using the incremental-iterative procedure mentioned 
previously. See references 1 and 3 for more detailed descriptions. 

All of the above steps and the required convergences are automatically per- 
formed in METCAN and are transparent to the user. 
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METCAN predicted longitudinal stress/strain responses for four ply 
[0/90/90/0] PlOO/Cu after a specific number of accumulated mechanical cycles 
are shown in figure 8. These curves were generated by assuming that: (1) the 

limit number of mechanical load cycles was 1 million, and (2) the cyclic stress 
amplitudes were the following ratios of the reference stress (static fracture 
stress): 85 percent for 100 000 cycles, 65 percent for 500 000 cycles, and 50 
percent for 800 000 cycles. It is interesting that the reference curve and the 
100 000 cycles curve exhibit some nonlinear behavior while the other two curves 
are practically linear. The important observations from the results shown in 
figure 8 are: (1) mechanical cyclic loading degrades the cyclic stress magni- 

tude but does not induce any discernable nonlinearities on the stress/strain 
behavior; and (2) the curves do not' start at zero strain because of residuals 
due to fabrication. Corresponding results for thermal load cycling are shown 
in figure 9. The number of accumulated thermal load cycles and the temperature 
range are shown in the figure. Note that the curves start at nonzero strains 
because of the fabrication residuals already mentioned and because the curve 
is generated at the elevated temperature. The limit number of thermal cycles 
was assumed to be 400 for generating these curves. Similarly METCAN predicted 
thermal cyclic load results are shown in figure 10 for a tungsten/copper HTMMC 
and in figure 11 for sill con-car bi de/ti tani um-al umi n i de - The important point 
to be observed from the thermal cycling load curves is that thermal cycling has 
severe degradation effects on the endurance of HTMMC 1 s. 


SIGNIFICANCE OF THE INTERPHASE AND IN SITU PROPERTIES 

HTMMC either have no metallurgical bond at the Interface or develop an 
Interphase as a result of the fiber-matrix reaction (interdiffusion) . The 
micromechanics model in METCAN incorporates the Interphase as a separate con- 
stituent with Independent properties. As such the effects of bond integrity 
on HTMMC behavior can be evaluated with judicious use of METCAN by using an 
Interphase with suitable properties. Results from one such evaluation are 
shown in figure 12. As would be expected, both the Interface and the in situ 
properties of the matrix have dramatic effects on the transverse tensile 
strength of this particular HTMMC. It is noted that the average experimen- 
tally measured strength is about 2 ks 1 which corresponds to the lowest value 
shown in figure 12. This value was obtained by assuming no Interface and 
annealed at yield in situ matrix conditions as inputs to METCAN. The important 
observation in figure 12 is that METCAN can be used to assess the fabrication 
process as well as the effects of variable constituent in situ behavior on 
HTMMC performance. For example, the lowest value indicates inadequate fabri- 
cation quality and inefficient use of copper as a matrix; while the highest 
value indicates that there is considerable room for improvement. This type of 
information is routinely obtained using METCAN, which provides an expedient 
capability for conducting a variety of parametric studies. 


SIMULATION OF COMBINED THERMOMECHANICAL CYCLING AND INTERFACE EFFECTS 

METCAN simulations of the effects of combined thermomechani cal cycling are 
shown in figure 13 for a SiC/Ti HTMMC. Comparisons with experimental data of 
the half cycle stress-strain response for a unidirectional composite and for an 
angleplied E+30D2S laminate are also shown. The dip in the experimental data 
represents, in part, the onset of excessive nonlinearities caused by high 
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Intralaminar shear stress in the 30° plies under displacement control testing. 
The experimental data were obtained by Dr. B. Lerch, Structures Division, NASA 
Lewis. The individual load cyclic effects are included for comparison pur- 
poses. As can be seen, combined cycling (lower right curves) degrades the 
mechanical behavior of HTMMC at a faster rate than do either of the individual 
loadings or their algebraic superposition. The important observation from the 
above discussion is that traditional isothermal cyclic data will be inadequate 
and even misleading for assessing the combined thermomechanical cyclic loading 
degradation effects on the behavior of HTMMC. 

The effects of an Interphase on thermal load cycling (fatigue), assuming 
400 maximum thermal cycles between 70 and 1000 °F, are shown in figure 14. 
Comparisons to those without the interphase at the end points are summarized 
in the fol lowing table. 


Number of Withou 

accumulated stres 

cycles ks 1 

1 78 

100 69 

200 69 

300 43 


t Interphase With 
s (strain), 
(percent) 

(0.44) 

(0.42) 

(0.42) 

(0.43) 


5-percent interphase 
stress (strain), 
ks i (percent) 

66 (0.44) 

60 (0.42) 

49 (0.42) 

40 (0.46) 


About the only observation to be made from these results is that the presence 
of the interphase tends to progressively degrade the stress magnitude at frac- 
ture, but not the fracture strain level for all practical purposes. Other 
tentative observations are: (1) the interphase has negligible effect on the 

fracture stress magnitude at relatively high thermal cycles, and (2) it has 
moderate degradation effects (about 16 percent) at relatively low number of 
thermal cycles. Again, the important conclusion is that METCAN can be used, 
rather routinely, to evaluate these effects which are complex and difficult to 
experimentally measure. Also METCAN can be used to aid the planning of stra- 
tegic experiments to selectively measure some of these effects. 


CONCLUDING REMARKS 

The significant conclusions of an Investigation to develop a procedure 
which can be used to computationally simulate the cyclic behavior of HTMMC are 
as follows: 

1. The cyclic behavior of HTMMC can be computationally simulated by using 
high temperature composite mechanics in conjunction with a multifactor- 
interaction relationship and with an incremental/interative nonlinear solution 
procedure. 

2. The computational simulation can be performed starting with the fabri- 
cation process, through thermomechanical cyclic loadings, and accounting for 
the Interphase. 

3. A computer code has been developed to perform these types of computa- 
tional simulations in an effective and routine manner. 
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4. Results from this computer code show that combined thermomechanical 
cycling degrades the HTMMC structural integrity at a faster rate than the 
algebraic superposition of the two Individual effects. 

5. The Interphase mainly affects the fracture stress at practically all 
levels of thermal cycling; it minimally effects the fracture strain, and then 
only at relatively high numbers of thermal cycles. 

6. The Interface has dramatic effects on the transverse tensile strength 
as do the in situ properties of the matrix. The ranges of values of the prim- 
itive variables, which contribute to these effects can be quantified. 

7. The type of computational simulation (or the respectively computer 
code) described herein can and should be validated and verified continuously 
in order to ascertain its predictive adequacy. 
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TABLE I. - ROOM TEMPERATURE CONSTITUENT PROPERTIES 
USED IN THE CURRENT STUDY 


Property type 

Fiber 


P100 

SiCA 

Borsic 

Tungsten 

T n a. 

6600 

4870 

6000 

6170 

E&i l > m P s ’ 

105 

62 

58 

59 

£2,22’ m P si 

0.9 

62 

58 

59 

G ftl2> m P si 

1.1 

23.8 

29 

22.7 

g 2.23> m P si 

0.7 

23.8 

29 

22.7 


0.2 

0.3 

0.2 

0.29 

v !l23’ in. /in. 

0.25 

0.3 

0.2 

0.29 

<*&]], ppm 

-0.9 

1.8 

2.8 

2.5 

a 422’ PP m 

5.6 

1.8 

2.8 

2.5 

Ka n , Btu/hr/i n . /°F 

7.5 

0.75 

1.9 

8.3 

K^22* Btu/hr/in./°F 

0.75 

0.75 

1.9 

8.3 

Sjm, Btu/hr/i n./°F 

325 

500 

510 

390 

Sfi. 1 1 C * Btu/hr/in. /°F 

200 

650 

600 

390 

s l22T» Btu/hr/in./°F 

25 

500 

100 

390 

S^22C* Btu/hr/in ./°F 

25 

650 

200 

390 

1 2S » Btu/hr/i n . /°F 

25 

300 

100 

236 

Sg^S* Btu/hr/i n ./°F 

12.5 

300 

100 

236 


Property type 

Matrix 


Copper 

Ti 

T °F 

‘mm* r 

1980 

1800 

E m , mpsi 

17.7 

12.3 

v m , in. /in. 

0.3 

0.32 

“m’ PP m 

9.8 

4.5 

K m , Btu/hr/in. /°F 

19.3 

0.39 

G mT ’ ksi 

32 

130 

S mC , ksi 

32 

130 

S mS , ksi 

19 

91 


















TABLE II. - GRAPHITE/COPPER ROOM TEMPERATURE MECHANICAL AND THERMAL PROPERTIES: 
COMPARISON OF METCAN PREDICTIONS AND 3-D FINITE ELEMENT ANALYSIS 
[PlOO/Cu at four different fiber volume ratios (FVR).] 


Property type 

FVR 

0.063 

0.2234 

0. 

466 

0.622 

METCAN 

3-DFEM 

METCAN 

3-DFEM 

METCAN 

3-DFEM 

METCAN 

3-DFEM 

E sm. m P s i 

23.7 

23.1 

37.3 

36.7 

58.4 

58.0 

72.0 

71.8 

Efi22’ m P si 

13.8 

14.9 

10.2 

11.0 

6.5 

6.9 

4.6 

5.0 

e 2.33> m P si 

13.8 

14.9 

10.2 

11.0 

6.5 

6.9 

4.6 

5.0 

Gq, 12 . mpsi 

5.8 

6.0 

4.5 

4.8 

3.2 

3.4 

2.5 

2.6 

Gjj,23> m P si 

5.6 

5.8 

4.2 

4.5 

2.8 

2.7 

2.1 

1.8 

G !ll3> m P si 

5.8 

5.8 

4.5 

4.5 

3.2 

2.7 

2.5 

1.8 

v 4l2* 

.29 

.30 

.28 

.30 

.25 

.29 

.24 

.28 

v^ 23 > in. /in. 

.30 

.30 

.30 

.26 

.30 

.22 

.30 

.19 

v !ll3» in. /in. 

.30 

.30 

.30 

.26 

.30 

.22 

.24 

.22 

PP m 

6.8 

6.6 

3.3 

3.7 

1 .0 

1 .40 

.18 

.42 

<*2,22’ PP m 

10.4 

10.5 

11.0 

10.9 

11.2 

11.1 

11.1 

10.9 

<*2,33’ ppm 

10.4 

10.5 

11.0 

10.9 

11.2 

11.1 

11.1 

10.9 

Kgj i , Btu-in./°F-hr-in.2 

19.7 

20.6 

20.6 

21.6 

22.0 

23.2 

22.8 

24.0 

K^22 j Btu-in./°F-hr-in.2 

15.6 

17.3 

11.8 

13.5 

7.8 

8.8 

5.8 

6.3 

<^23* Btu-in./°F-hr-in.2 

15.6 

17.3 

11.8 

13.5 

7.8 

8.8 

5.8 

6.3 


SYSTEM 

EQUATIONS 

OF 

MOTION 


[M] (ii) + [C] (u] + [K](u) = (F(t)] 


^ (■) * (Ua) 


CONSTITUTIVE 

EQUATION 


M = [D][B]fuj 


V (cl < (SaJ 


NATURAL 

FREQUENCY 

EIGENPROBLEM 


([K]-u, 2 [M])[uS = [o] 


U) < b>A 


FIGURE 1. - OVERVIEW OF STRUCTURAL ANALYSIS. 
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FIGURE 2. - METAL-MATRIX COMPOSITE FABRICATION PROCESS. 
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FIGURE 3. - INTEGRATED MULTI-SCALE COMPUTATIONAL SIMULATION. 
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RATIONALE: 


SUBREGIONS OF 

INTRALAMINAR 

NONUNIFORMITY 



r- d-H 


• GRADUAL EFFECTS DURING MOST RANGE, RAPIDILY DEGRADING NEAR FINAL STAGES 

• REPRESENTATIVE OF THE IN SITU BEHAVIOR FOR FIBER, MATRIX, INTERPHASE, COATING 

• INTRODUCTION OF PRIMITIVE VARIABLES (PV) 

• CONSISTENT IN SITU REPRESENTATION OF ALL CONSTITUENT PROPERTIES IN TERMS OF PV 

• ROOM TEMPERATURE VALUES FOR REFERENCE PROPERTIES 

• CONTINUOUS INTERPHASE GROWTH 

• SIMULTANEOUS INTERACTION OF ALL PRIMITIVE VARIABLES 

• ADAPTABILITY TO NEW MATERIALS 

• AMENABLE TO VERIFICATION INCLUSIVE OF ALL PROPERTIES 

• READILY ADAPTABLE TO INCREMENTAL COMPUTATIONAL SIMULATION 


NOTATIONS: 

P - PROPERTY; T - TEMPERATURE; S - STRENGTH; R - METALLURGICAL REACTION; N - NUMBER OF CYCLES; t - TIME; 
OVER DOT - RATE; SUBSCRIPTS: 0 - REFERENCE; F - FINAL; M - MECHANICAL; T - THERMAL 

FIGURE 4. - MULTIFACTOR INTERACTION RELATIONSHIP FOR IN SITU CONSTITUENT MATERIAL BEHAVIOR. 




FIGURE 5. - UNIDIRECTIONAL TUNGSTEN/COPPER MMC ROOM TEM- 
PERATURE STRESS-STRAIN RESPONSE. 


FIGURE 6. - Gr/Cu COMPOSITES ROOM TEMPERATURE STRESS- 
STRAIN BEHAVIOR (UNIDIRECTIONAL LAMINATE AT .4 FVR). 
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STRESS EPS XX, PERCENT 


FIGURE 7. - UNIDIRECTIONAL BORS I C/T ITANIUM MMC ELASTIC 
MODULI AT ROOM TEMPERATURE. 


FIGURE 8. - MECHANICAL CYCLIC LOAD EFFECTS (PIOO/Cu AT .4 
FVR) . 



STRAIN EPS XX, PERCENT 


FIGURE 9. - THERMAL CYCLIC EFFECTS OF PIOO/Cu COMPOSITE 
(AT .4 FVR). 



STRAIN EPS XX, PERCENT 

FIGURE 10. - THERMAL CYCLIC EFFECTS OF TUNG/Cu COMPOSITE 
(AT .4 FVR). 
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FIGURE 11. - THERMAL CYCLIC EFFECTS OF SiCa/Ti 15 COMPOSITE 
(AT .H FVR). 



YIELD FRACTURE YIELD FRACTURE 

FIGURE 12. - UNIDIRECTIONAL GRAPHITE/COPPER MMC TRANSVERSE 

STRENGTH BOUNDS. 
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